f=[3,12,3,4,0,11,2,5,9,0,6,7,1,5,0]';
%intcon=[1 2 3 4 5 6 7 8 9 10 11 12 13 14 15];
A=[];b=[];
Aeq=[1,1,1,1,1,0,0,0,0,0,0,0,0,0,0;
     0,0,0,0,0,1,1,1,1,1,0,0,0,0,0;
     0,0,0,0,0,0,0,0,0,0,1,1,1,1,1;
     1,0,0,0,0,1,0,0,0,0,1,0,0,0,0;
     0,1,0,0,0,0,1,0,0,0,0,1,0,0,0;
     0,0,1,0,0,0,0,1,0,0,0,0,1,0,0;
     0,0,0,1,0,0,0,0,1,0,0,0,0,1,0;
     0,0,0,0,1,0,0,0,0,1,0,0,0,0,1];

    beq=[8;5;9;4;3;5;6;4];
    lb=[0 0 0 0 0 0 0 0 0 0 0 0 0 0 0];
    [x,yunfei]=linprog(f,A,b,Aeq,beq,lb);
    x=reshape(x,5,3);
    x=x'
    yunfei
    
    
    

